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FOREWORD 


The  Army  Multiple  Systems  Evaluation  Program  (MSEP)  is  a 
comprehensive  program  developing  general  analytic  techniques  for  the 
prediction  of  high  electromagnetic-pulse  vulnerability  and  heurdening 
technology  and  for  the  application  of  these  techniques  to  a list  of 
critical  systems.  The  2uialytic  techniques  have  been  verified  for  a 
large  class  of  tactical  systems.  The  hardening  techniques  have  been 
applied  to  specific  systems  and  are  now  resulting  in  product  improvement 
programs  leading  to  hardened  equipment  in  the  field. 
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1 . INTRODUCTION 

This  effort  was  written  under  the  sponsorship  of  the  Multiple  Systems 
Evaluation  Program  (MSEP) , which  has  as  its  main  objective  to  harden 
Army  tactical  systems  to  the  electromagnetic  pulse  (EMP)  associated  with 
the  exoatmospheric  detonation  of  a nuclear  weapon.  Along  with  this 
major  objective,  MSEP  is  tasked  also  with  the  aim  to  develop 
experimental  and  analytical  evaluation  techniques  that  are  applicable  to 
all  systems'  problems.  These  objectives  are  satisfied,  since  SEMCON 
provides  an  analytic  method  to  examine  various  semiconductor  devices 
used  in  Army  systems  and  assists  in  making  recommendations  for  hardening 
modifications. 

When  making  EMP  vulnerability  assessments,  it  is  often  necessary  to 
perform  a circuit  analysis  on  specific  circuits  driven  by  a voltage 
induced  by  EMP.  This  analysis  often  points  out  a large  power  value 
produced  on  various  semiconductor  devices  that  could  cause  the  devices 
to  fail  if  sufficient  power  were  present.  However,  to  know  if  a 
particular  device  fails  and  at  what  power  level,  a semiconductor  device 
must  be  analyzed  and  a damage  curve  fitted  to  certain  data.  This  curve 
gives  a good  indication  as  to  the  vulnerability  of  the  device  to 
particular  thresholds  of  power.  Therefore,  the  main  purpose  of  computer 
code  SEMCON  is  to  provide  the  necessary  information  about  semiconductor 
devices  so  that  it  is  possible  to  make  various  EMP  vulnerability 
assessments . 

The  semiconductor  damage-data-reduction  computer  code,  SEMCON, 
analyzes  various  semiconductor  devices.  The  code  was  written  for  a 
Control  Data  Corporation  (CDC)  6600  series  computer  system  located  at 
the  Mobility  Equipment  Research  and  Development  Center  (MERDC) , Fort 
Belvoir,  VA.  Written  to  run  employing  the  FTN  compiler,  SEMCON  is 
operational  using  the  SCOPE  3.4.3  control  language.  To  evaluate  a 
semiconductor  device,  SEMCON  is  capable  of  giving  the  following 
information : 

a.  Determination  of  damage — the  failure  or  nonfailure  of  a device 

b.  Calculations  of  power,  energy,  and  impedance 

c.  Least-squares  fits  of 
p = at“ ^ , 

P = bt"^, 
p = at~^  + bt~*5 

to  either  power  versus  failure  time  or  power  versus  pulse  width 
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d.  Selection  of  the  best  least-squares  fit  to  specified  data 

e.  Log-log  plots  of  least-squares  fits  in  either  the  forward  or 

reverse  direction  ] 

f.  Log-log  plots  of 

TD  versus  Z — failure  time  versus  impedance  ' 

TD  versus  VOC — failure  time  versus  open-circuit  voltage  • 

TDOC  versus  VOC — pulse  duration  versus  open-circuit  voltage 

« J 

TD  versus  IP — failure  time  versus  peak  current  j 

in  either  the  forward  or  reverse  direction  j 

g.  Tabular  listing  of  input  data  and  all  calculations  | 

i 

The  subsequent  sections  detail  the  development  and  capabilities  of  ; 

SEMCON,  along  with  the  necessary  information  to  utilize  the  code. 

i 

2.  FAILURE  DETERMINATION  AND  POWER,  ENERGY,  AND  IMPEDANCE  CALCULATIONS  I 

To  determine  the  damage  on  a particular  pulse  number,  the  specific 
values  used  in  the  calculations  are  (1)  VZBP  and  VZAP,  (2)  VZBN  and 
VZAN,  and  (3)  GAINB  and  GAINA.  The  criteria  used  to  determine  failure 
are  any  one  of  the  following  three  conditions: 

a.  .8*VZBP  - VZAP  > 0 

b.  .8*VZBN  - VZAN  > 0 

C.  .8*GAINB  - GAINA  >0 

If  any  one  of  these  conditions  is  satisfied,  then  for  that  pulse  number  ' 

the  device  is  said  bo  fail. 

Another  of  the  features  of  SEMCON  is  the  calculations  of  the  power 
(P) , energy  (E) , and  impedance  (Z) , These  computations  are  very 
straightforward  and  are  found  by 

P = VP*TD, 

E = P*TD, 

VP  - VZBP 
^ ■ — 1? — • 
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3.  DERIVATION  OF  FITTING  PARAMETERS 

The  main  purpose  of  SEMCON  is  to  fit  various  curves  to  specified 
data.  They  are  fit  by  using  the  method  of  least  squares  to  calculate 
the  various  fitting  parameters.  The  plots  of  these  curves,  as  well  as 
the  data,  are  scaled  logarithmically,  since  there  can  be  a wide  range  in 
the  values  of  the  data.  Thus,  it  is  necessary  to  calculate  the 
least-squares  fits  in  terms  of  logarithms,  which  introduces  a nonlinear 
problem.  Two  methods  are  used  to  calculate  the  least-squares  fits  in 
SEMCON.  One  uses  the  minimizing  conditions  known  as  the  normal 
equations,  and  the  other  uses  a minimization  technique  to  find  the 
minimum  of  a function  of  two  variables.  In  the  following  paragraphs, 
the  two  methods  are  explained  fully  along  with  a detailed  derivation  of 
the  least-squares  fitting  parameters. 

The  first  method  used  to  find  the  least-squares  fits  involves  the 
solution  of  the  normal  equations.  This  technique  involves  substituting 
the  n data  pairs  into  the  fitting  equation. 


f (x)  = ai(|>i  (x)  + a2i(>2(x)  + . . . + 


where  <j>.  (x)  are  known  functions  of  x,  and  a.  are  the  k linear  fitting 
parameters,  thus  obtaining  n linear  equations  with  k unknowns.  The  next 
step  is  to  multiply  each  equation  by  its  corresponding  coefficient  of 
the  first  unknown,  which  results  in  n new  equations.  Adding  these 
equations  gives  the  first  normal  equation.  This  process  is  repeated  k 
times  with  respect  to  each  of  the  unknown  parameters,  so  that  we  end  up 
with  k equations  in  k unknowns  (the  fitting  parameters) . The  resulting 
system  of  simultaneous  linear  equations,  known  as  the  normal  equations, 
is  then  solved  for  the  k fitting  parameters. 

In  our  particular  case,  we  wish  to  use  this  method  to  fit 

p = at“^, 

p = bt“^ 

to  experimental  data,  where  p is  the  power,  t is  the  failure  time,  and  a 
and  b are  the  fitting  parameters.  Since  plots  are  desired  that  are 
scaled  logarithmically,  we  must  find  the  parameters  a and  b by  applying 
the  least-squares  technique  to 


log  p = log  af^. 

(1) 

log  p = log  bt“\ 

(2) 
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All  logarithms  are  assumed  to  be  to  the  base  10.  By  the  elementary 
properties  of  logarithms,  we  clearly  have  a linear  problem,  since 

log  p = log  a - log  t, 

log  p = log  b - *5  log  t. 

First,  we  consider  log  p = log  at~^  and  suppose  we  have  n data  pairs 
(ti,  p^)  for  1 - i - n.  Substituting  into  equation  (1) , we  get 

log  pj  = log  at“^, 

log  p^  = log  at~  ^ , 


log  p = log  at  ^ , 
n n 

and  simplifying  gives 

log  Pj  = log  a - log  tj , 

log  ?2  = log  a - log  t^. 


log  p = log  a - log  t . 
n n 

Since  this  system  is  linear  in  terms  of  logarithms,  multiplying  by  the 
coefficient  of  log  a in  each  equation  and  adding  gives 


n n n 

Yi  ^05  Pi  = Y log  a - Y 1°9  t . 
i=l  i=l  i=l 


This  is  the  only  normal  equation,  since  there  is  only  one  fitting 
parameter.  Solving  for  log  a,  we  obtain 


n n 

Y Pi  + Z!  i°5  ^i 

i=l  i=l 


log  a = 


n 


i 


which  gives  the  solution  for  a as 


a = 10  , 


where 


X = 


1 / n n \ 

^ ""ij' 


Now,  considering  log  p = log  bt  and  proceeding  as  above,  we  get 
log  Pi  = log  b - *j  log  t^  , 
log  p^  = log  b - *3  log  t^  , 


log  p = log  b - *j  log  t 
n n 


This  system  also  is  linear  in  terms  of  logarithms,  so  multiplying  by  the 
coefficient  of  log  b in  each  equation  and  adding  gives 

n n n 

^ log  P.  = log  b - >5  ^ log  t . 
i=l  i=l  i=l 


Again,  this  is  the  only  normal  equation,  since  there  is  only  one 
fitting  parameter.  Solving  for  log  b,  we  obtain 


log  b = 


n n 

log  p.  + *5  X]  t 

i=l  ^ i=l  ^ 


which  gives  the  solution  for  b as 


where 


b = 10^, 


1 /n  n ) 

Y = -(  Z + 4 L t 

V=i  J 
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The  second  method  employed  to  calculate  the  least-squares  fits 
involves  finding  a minimum  of  a function.  Suppose  we  want  to  fit 

f (x)  = ai4>i(x)  + a24^2<^) 

to  n data  points,  where  ♦^(x),  <t>2  (x)  are  know  functions  of  x,  and  , 
are  the  fitting  parameters.  Let  fj^  = f(x.)  for  1 i i 5 n be  the 

experimental  data  that  we  wish  to  fit.  Then  by  the  principle  of  least 
squares,  we  want  to  minimize 

F(ai,  a2>  = X!  w(x^)|^f^  - , 

where  w(x^)  is  the  weight  function  for  each  x^.  For  our  case,  we  choose 

w (x , ) = 1 for  1 5 i £ n . 

1 

Thus,  we  must  find  the  minimum  of 


F(a, 


^2^ 


+ 


To  find  the  minimum,  we  proceed  in  the  following  manner: 


a.  Start  with  initial  values  a j , a^ 

b.  At  the  ith  step,  halve  aj^^^  , and  try  to  step  first  in  the 
positive  and  then  in  the  negative  x-direction,  to  see  if  either 
resulting  value  gives  a minimum,  i.e. , 


^ (i) 

(i)  ^ (i) 

(i) 

(i) 

«■  F 

ai 

, 32 

c.  Halve  a2^  , and  try  to  step  first  in  the  positive  and  then  in 
the  negative  y-direction,  to  see  if  either  resulting  value 
gives  a minimum,  i.e.. 


F aj  , 32  I 


<F  4^^ 


If  neither  (b)  nor  (c)  gives  a minimum  halve  a{^^  and  a^^^ 
again,  and  repeat  (b)  and  (c) , i.e.,  j = 2. 

If  a minimum  results  from  either  (b)  or  (c) , then  let 


(i+1)  (i)  ^ (i+1)  (i)  . , , 

Lj  = aj  ± r-  #32  = 32  , if  only  (b)  holds, 


44^ 

= aj^^  ± ^ = a|^^  , if  only  (c)  holds, 


(i)  ^ ^(i+1)  _ ^(i)  ^ . 

= aj  ± r I ^2  ~ ^2  - T ' 

oJ 


if  (b)  and  (c)  hold. 


and  check 


(i+1)  (i+1)]  (i)  (i)] 

F|^ai  ,32  J - Fj^ai  , 32  J 


where  € is  a chosen  convergence  tolerance.  If  equation  (3)  holds,  then 


, ^ (i+1)  (i+1) 

aj  + a]^  f a2  “ ^2 


are  the  fitting  parameters.  If  equation  (3)  does  not  hold,  the  process 
is  repeated. 

This  technique  is  very  simple  and  always  leads  to  a minimum 
although  it  may  l)e  only  a local  minimum.  If  a local  minimum  is  reached 
which  is  not  close  to  the  true  minimum,  then  it  seems  reasonable  to 
assume  that  the  least-squares  fit  will  be  poor.  Thus,  when  the  fit  is 
bad,  it  could  probably  be  attributed  to  this  distant  local  minimum, 
although  there  are  other  occurrences  that  could  cause  a bad  fit.  The 
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only  solution  to  this  problem  is  to  try  the  minimization  technique  again 
with  different  initial  values.  A repetition  of  this  process  leads  to 
the  true  minimum  of  the  function  F(aj,  32).  The  above  technique  was 
limited  to  two  fitting  parameters,  but  it  is  an  easy  matter  to  extend 
this  method  to  k parameters. 

For  our  purpose,  we  wish  to  use  this  method  to  fit 


-1 

p = at  + bt  (4) 

to  experimental  data  (t.,  p.)  for  1 S i < n,  where  p is  the  power,  t is 
the  failure  time,  n is  the  number  of  data  points,  and  a and  b are  the 
fitting  parameters.  This  poses  a nonlinear  problem,  since  log 
(at“^  + bt~**)  cannot  be  linearized  as  was  log  (at~M  and  log  (bt~^). 
This  problem  is  easily  overcome  by  using  the  minimization  technique 
described  above.  Let 


P(t;  a,b)  = log  (at“^  + bt”*’*)  , 
and  we  want  to  minimize 


S(a,b)  = ^ w(t^)  jlog  p^  - > 

where  w(t.)  is  the  weight  function  for  each  t. . Since  we  choose 
X 1 


w(t^)  = 1 for  i = 1,  2,  . . . n. 


we  must  minimize 


S(a,b)  = ^ |log  p^  - log  (at?  ^ 

The  minimization  for  S(a,b)  is  identical  to  the  minimization  technique 
given  above.  A flow  chart  for  finding  the  minimum  of  S(a,b)  and  thus 
the  fitting  parameters  a and  b is  given  in  figure  1. 


Although  the  techniques  outlined  above  are  simple  and 
straightforward,  caution  should  be  taken  when  interpreting  the 
least-squares  fits.  First,  the  fits  are  only  as  good  as  the  data. 
Poorly  taken  and  recorded  data  result  in  meaningless  fits.  Probably  the 
major  fault  with  least  squares  is  that  a single  very  wrong  measurement 
greatly  distorts  the  results,  because  in  the  squaring  process,  large 


a.b  OBTAMEO  FROM 


L5 


•(q'B)s  BuTZTmTUTUi  JOJ  *T  sjnbxj 


residuals  are  dominant — one  gross  error  10  times  larger  than  most  of  the 
others  has  the  same  effect  in  the  sum  of  the  squares  as  have  100  of  the 
others.^  Second,  a considerable  amount  of  mathematical  and  statistical 
sophistication  can  be  employed  to  give  the  fits  much  more  reliability 
when  the  proper  data  are  involved.  Several  statistical  techniques  can 
be  applied  to  these  least-squares  fits  that  would  enhance  their 
assurance.  One  is  a statistical  test  of  hypothesis  concerning  the 
confidence  in  the  fitting  parameters  a,  b in  equations  (1) , (2) , and  (4) . 
The  other  would  be  to  test  the  goodness  of  fit  by  using 
the  statistical  test.  These  two  techniques  can  be  utilized  to  make 
SEMCX»N's  least-squares  fits  more  reliable. 

4.  GAUSS  CRITERION  OF  GOODNESS  OF  FIT 


One  of  the  main  questions  pertaining  to  least-squares  fits  concerns 
the  goodness  of  fit  of  the  curve  to  the  data.  This  problem  can  be  dealt 
with  when  using  the  least-squares  technique,  cind  probably  the  most 
widely  used  method  is  the  X goodness  of  fit  statistical  test.  Although 
the  X^  test  has  not  been  utilized  in  SEMCON,  a test  to  determine  which 
of  the  curves 


p = bt~** 
p = afl  + bt“*5 


best  fits  the  data  is  applied.  This  test,  toown  as  the  Gauss  criterion 
of  goodness  of  fit,  is  now  explained.  Let  p^^  be  an  experimental  value, 
p.  a value  computed  from  a least-squares  fit,  n the  number  of  data 
points,  and  m the  number  of  fitting  parameters  in  the  relation.  Then  we 
define 


n 


and  the  Gauss  criterion  of  goodness  of  fit  states  that  the  best  fit  is 
the  one  that  minimizes  Q.  Since  the  plots  of  these  fits  are  scaled 
logarithmically,  we  calculate 


^R.  W.  Hamming,  Numerical  Methods  for  Scientists  and  Engineers , 2nd 
edition,  McGraw-Hill  Book  Co.,  Inc.,  New  York  (1973). 
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and  find  the  curve  that  minimizes  Values  of  (i  are  printed  in  the 

output  of  SEMCON,  as  well  as  the  curve  that  best  fits  the  data.  This 
can  be  seen  in  appendices  A euid  B. 

5.  PLOTTING  INFORMATION 

There  are  several  ways  to  obtain  the  plots  generated  by  SEMCON. 
Through  MERDC,  SEMCON  has  the  capability  to  use  a Cal  Comp  model  835 
cathode  ray  tube  (CRT)  plotting  system.  This  system  takes  plotting 
information  from  a magnetic  tape  and  converts  the  data  to  be  plotted  to 
incremental  plotter  conmands,  which  produce  discrete,  electron-beam 
deflections  (relative  to  the  x and  y axes)  and  intensity  variations  on 
the  face  of  a CRT.  The  CRT  display  is  transmitted  through  a camera  lens 
system  and  recorded  on  35-mm  microfilm.  The  exposed  film  is  then 

processed,  and  the  plots  are  recorded  on  a magnetic  tape  at  the  MERDC 
computer  center.  A call  to  the  MERDC  computer  center  is  necessary  to 
initiate  the  hard-copy  processing  of  the  plots  from  the  magnetic  tape. 

In  addition  to  the  CRT  plotting  system,  SEMCON  has  the  ability  to 
use  a Mohawk  Data  Sciences  (MDS)  model  2400  remote  batch  terminal 
accompanied  with  a Houston  Instruments  COMPLOT  plotting  system.  This 
system,  located  at  the  Harry  Diamond  Laboratories  (HDL)  Woodbridge 
Research  Facility  (WRF) , enables  all  plots  and  output  information  to  be 
received  on  site. 

Each  of  the  methods  to  obtain  plots  from  SEMCON  has  different 
plotting  software.  Examples  of  this  software  are  contained  in 
appendices  A and  B. 

6.  OPTIONS  OF  SEMCON 

The  output  of  SEMCON  can  be  varied  substantially,  so  an  analysis  of 
the  possible  options  for  outputing  various  information  is  essential; 
The  following  is  a description  of  the  options  available. 

a.  Tabular  data  options 

(1)  Listing  of  tabular  data  only — no  plots  are  given. 

(2)  Listings  of  tabular  data  and  plots  are  given. 

(3)  Plots  are  given,  but  the  tabular  listing  is  suppressed. 
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Plotting  selections 


(1)  Available  items  versus  item  plots  are 

(a)  TD  versus  Z — failure  time  versus  impedance 

(b)  TD  versus  VOC — failure  time  versus  open-circuit 
voltage 

(c)  TDOC  versus  VOC — pulse  duration  versus  open-circuit 
voltage 

(d)  TD  versus  IP— failure  time  versus  peak  current. 

(2)  Plots  in  the  reverse  direction  only 

(3)  Plots  in  the  forward  direction  only 

(4)  Plots  in  both  forward  and  reverse  directions 
(.5)  No  item  versus  item  plots 

(6)  Reverse  direction  item  versus  item  plots  only 

(7)  Forward  direction  item  versus  item  plots  only 

(8)  Forward  and  reverse  direction  item  versus  item  plots  only 

(9)  Forward  and  reverse  direction  item  versus  item  plots  and 
least-squares  fits 

(10)  Plots  may  be  received  on  the  MDS  remote  batch  terminal  or 
at  MERDC. 

(11)  Plotting  software  may  be  either  the  regular  Cal  Comp 
plotting  software  or  the  software  available  through  the  MDS 
remote  batch  terminal. 

Input  data — The  regular  input  data  may  be  entered  or  a simpler 
form  of  input  data  that  requires  only  pulse  width  versus  power 
entires. 


Least-squares  fits — Three  curves  may  be  fitted  to  either 
failure  time  (TD)  versus  power  (P)  or  pulse  width  versus  power. 

(1)  Fit  to  p * at~^  and  p = bt~**  only. 
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(2)  Fit  to  p = at  ^ + bt  **  only. 

(3)  Fit  to  all  three  equations. 

7.  DESCRIPTION  OF  SEMCON  AND  ITS  SUBROUTINES 

The  SEMCON  code  consists  of  a main  program  and  19  subroutines.  The 
documentation  of  SEMCON  is  detailed  by  the  following  analysis  of  the 
main  program  and  its  subroutines.  A listing  of  SEMCON  can  be  found  in 
appendix  C. 

SEMCON  reads  in  the  data  and  writes  the  tabular  output  seen  in 
appendix  A.  It  also  makes  the  power,  energy,  and  impedance 
calculations. 

FAILUR  calculates  the  failure  or  nonfailure  of  each  pulse  number. 

FRFAIL  separates  the  data  with  respect  to  failure  in  the  reverse 
or  forward  directions. 


DIRECT  sets  up  the  variables  T and  PP,  which  are  used  in  the 
4 calculations  of  the  three  least-squares  fits,  with  respect 

to  either  the  forward  or  reverse  direction.  Also,  DIRECT 
chooses  which  least-squares  fit  will  be  calculated. 

WRITE  outputs  the  various  fitting  parameters  for  the  chosen 
least-squares  fit.  Also,  some  other  useful  information  is 


written,  such  as  the  curve  that  best  fits  the  data  and  the 
minimum  of  the  function  S(a,b)  (sect.  4). 

PLOTT  contains  all  the  plot  titles  and  axes  information  and 
labels.  This  subroutine  initiates  the  plotting  of  all 
curves  by  calling  PWRITE. 

ALSQAR  calculates  the  fitting  parameter  for  p = at  ^ . 

BI£QAR  calculates  the  fitting  parameter  for  p = bt"*'’  . 

ABLSQAR  calculates  the  fitting  parameters  for  p = at  ^ + bt  . 

PWRITE  performs  the  plotting  using  the  Cal  Comp  model  835  CRT 
plotting  system  by  calling  the  WRF  library  routines  in 
TALPLOT.  The  subroutines  SELPLOT,  TALGRAF,  and  TALDATA  of 
TALPLOT  are  referenced  in  PWRITE  only. 

SELPLOT  selects  the  beam  intensity  of  the  Cal  Comp  model  835 
microfilm  plotter  and  also  enables  the  user  to  choose  if  he 
wants  the  plots  sent  to  the  MDS  system  at  WRF. 
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TALGRAF  sets  up  the  plot  axes  and  labels  and  specifies  the  number 
of  sets  of  data  to  be  plotted  on  one  graph. 

TALDATA  plots  a set  of  (x,y)  data  pairs. 

SORT  sorts  the  data  in  array  T of  subroutine  DIRECT  in 

increasing  order  so  that  the  plotting  is  done  correctly. 

GOODFIT  calculates  the  sum  of  the  squares  of  the  difference  between 
the  experimental  value  and  the  corresponding  value  computed 
from  a least-squares  fit  for  each  least-squares  fit.  This 
sum  is  used  in  the  Gauss  criterion  of  goodness  of  fit. 


GAUSS  Determines  the  best  fit  for  the  data  by  the  Gauss  criterion 
and  outputs  the  correct  fit. 

SAB  calculates  the  functional  values  of  S{a,b)  (sect.  4). 


OWNDAT  enables  the  user  to  enter  only  data  for  pulse  width  versus 
power  input.  It  outputs  these  data  by  calling  COLMNS  eind 
makes  the  appropriate  calls  to  the  other  routines  to 
calculate  the  least-squares  fits  and  make  the  desired 
plots. 

COLMNS  outputs  the  data  for  pulse  width  versus  power  in  increasing 
columns. 


CHANGE  changes  the  abscissa  and  ordinate  titles  to  correspond  with 
the  use  of  data  for  pulse  width  versus  power.  It  also 
enters  those  data  into  the  plotting  software. 


8.  DATA  INPUT  PREPARATION 


Input  cards  are  prepared  for  SEMCON  in  the  following  manner. 


Card  Columns  Variable  name  Format 


Comment 


Card  set  1:  Plot  and  parameter  card 

1 10  IND  II 


IND  = 1,  plot  in  reverse 
direction  only 
IND  = 2,  plot  in  forward 
direction  only 
IND  = 3,  plot  in  both 
directions 
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Card 

1 


Columns 
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Variable  name 


FLAG 


Format  Comment 


II  FLAG  = 1,  least-squares 
fit  to  p = at~^  and 
p = bt“T 

FLAG  = 2,  least-swares  fit 
to  p = at”  + bt”^ 

FLAG  = 3,  least-squares  fit 
to  all  three  equations 


30  IPRAM  II  IPRAM  = 0,  no  item  versus 

item  plots 

IPRAM  = 1,  reverse  direction 
item  versus  item  plots  only 
IPRAM  = 2,  forward  direction 
item  versus  item  plots  only 
IPRAM  = 3,  forward  and 
reverse  direction  item 
versus  item  plots  only 
IPRAM  = 4,  forward  and 
reverse  direction  item 
versus  item  plots  and  all 
curves  as  specified  by 
IND  and  FLAG 


40  ISTOP 


49-50  IPLOT 


60  ISEM 


11  ISTOP  = 1,  tabular  output  only 
ISTOP  = 2,  tabular  output  and 

all  output  designated  by  IND, 
FLAG,  and  IPRAM 
ISTOP  = 3,  no  tabular  output, 
but  all  plots  as  indicated 

12  IPLOT  = 0,  plots  come  out 

on  MDS  system 

IPLOT  = K,  vdiere  1 < K < 36 
indicates  beam  intensity 
of  Cal  Cc»np  plotter;  a good 
value  is  IPLOT  = 18;  plots 
come  out  at  MERDC 
NOTE:  ITYPE  must  be  equal  to  0. 

II  ISEM  = 0,  enter  all  input 
data  in  card  sets  3 and  4 
ISEM  = 1,  enter  only  pulse 
width  versus  power  data  on 
card  set  5 
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Card 

Columns 

Variable  name 

Format 

Comment 

1 

70 

ITYPE 

II 

ITYPE  = 0,  regular  Cal  Comp 
plotting  software  is  used 
ITYPE  = 1,  plotting  software 
available  on  MDS  remote 
batch  terminal  is  used 

Card  set 

2;  Title 

card 

2 

1-3 

NUMB 

13 

If  ISEM  = 0 on  card  1, 

NUMB  is  total  number  of 
devices 

If  ISEM  = 1 on  card  1, 

NUMB  is  number  of  pulse 
width  versus  power  points 

4-11 

DEV 

A8 

Device  riome 

12-31 

DEVTYPl , DEVTYP2 

2A10 

Device  type 

32-40 

JUNC 

A9 

Junction  of  device 

41-60 

MANl,  MAN2 

2A10 

Manufacturer  of  device 

61-70 

TECH 

AlO 

Technician 

71-80 

DATEE 

AlO 

Date 

NOTE:  If  ISEM  = 1 on  card  1, 
skip  to  card  set  5. 

Card  set 

3;  Pulse 

card 

3 

1-2 

NPULS 

12 

Total  number  of  pulses 
for  particular  device 
number 

3-80 

- 

- 

Any  comments  or  notations 

may  be  placed  in  these  columns 
NOTE;  Card  set  3 is  repeated 

1 as  many  times  as  indicated  on 

card  set  2 for  variable  NUMB, 
and  thus  card  set  4 is  repeated 

I this  same  number  of  times. 
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Card  Columns 

Variable  neune 

Format 

Comment 

Card  set  4:  Parameter  value  cards 

4 1~3 

NDEV 

13 

Device  number:  total  number 
of  repetitions  of  this  card 
must  not  exceed  600 

4~5 

PULSE 

12 

Pulse  number;  not  to  exceed  ' 

6-15 

TDOC 

E10.3 

Pulse  duration 

16-25 

VOC 

F10.3 

Open-circuit  voltage 

26-35 

VP 

F10.3 

Peak  voltage 

36-45 

IP 

F10.3 

Peak  current 

46-55 

TD 

E10.3 

Failure  time  (TD  0) 

56-65 

VZBP 

F10.3 

Voltage  before  forward 
breakdown  voltage 

5 1-10 

VZAP 

P10.3 

Voltage  after  forward 
breakdown  voltage 

11-20 

VZBN 

F10.3 

Zener  voltage  before 
reverse  breakdown  voltage 

21-30 

VZAN 

F10.3 

Zener  voltage  after 
reverse  breakdown  voltage 

31-40 

GAINB 

ElO.  3 

Current  gain  before  device 
was  pulsed 

41-50 

GAINA 

E10.3 

Current  gain  after  device 
was  pulsed 

•“•51-52 

PIN 

A2 

Polarity  of  applied  pulse: 
A+  = forward  direction 
C+  = reverse  direction 

53-71 

REMAR1.REMAR2 

A10,A9 

Remarks 

NOTE:  Card  set  4 represents 
data  for  one  pulse  and  is 
repeated  as  many  times  as 
indicated  on  card  set  3 for 
variable  NPULS. 

If  ISEM  = 0 on  Card  1,  data 
input  is  complete 
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Card  Columns 

Variable 

name 

Format 

Comment 

Card  set  5:  Pulse 

width  and 

power 

cards 

6 1-20 

IDENT1,IDENT2 

2A10 

Identification 

21-30 

- 

- 

Blank 

31-40 

XLAB 

AlO 

Abscissa  label 

41-50 

- 

- 

Blank 

51-60 

YLAB 

AlO 

Ordinate  label 

7 1-10 

T 

E10.3 

Pulse  width 

11-20 

PP 

E10.3 

Power 

NOTE:  Data  for  pulse  width 
versus  power  must  be 
arranged  in  increasing 
order  with  respect  to 
pulse  width. 

Card  7 is  repeated  according 
to  value  NUMB  appearing  on 
Card  2 


A sample  listing  of  the  input  cards  when  ISEM  = 0 on  card  set  1 is 
given  in  appendix  D.  Also,  a sample  listing  of  the  input  cards  when 
ISEM  = 1 on  card  set  1 is  contained  in  appendix  E. 


9.  CONTROL  CARDS  FOR  RUNNING  SEMCON 

The  SEMCON  code  is  operational  using  the  SCOPE  3.4.3  control 
language  on  the  CDC  6600  computer  at  MERDC.  Since  several  different 
plotting  packages  can  be  used,  three  distinct  control  card  sets  are 
required  for  the  various  options  available.  The  following  details  the 
operation  of  the  code. 

If  IPLOT  = K where  IS  K S 36  on  card  set  1 (plots  come  out  at  MERDC 
using  the  Cal  Comp  plotting  software) , then  the  control  cards  for 
running  SEMCON  are  as  follows: 

EM (MTl,  T300) 

TASK(TNEM ,PW ,TRTS)  [user's  name] 

REQUEST,  TAPE50,  HI,  VSN=  ,RING. 
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NOTE:  VSN=  is  the  number  of  the  blank  tape  at  MERDC. 

ATTACH , AGO , B INS EMCON , ID=EM7 1 606 . 

ATTACH , F 1 , TALWY ATT ,ID=EM71602, MR= 1 . 

LIBRARY (Fl) 

MAP (PART) 

LDSET (PRESET=NGINF) 

AGO. 

TALP  LOT  (,  taped 
EXIT. 

7/8/9 

[Data] 

0/6/7/8/9 


j If  IPLOT  = 0 and  ITYPE  = 0 on  card  set  1 (plots  come  out  on  the  MDS 

system  using  the  Cal  Comp  plotting  software) , then  the  control  cards  for 
running  SEMCON  are  as  follows: 

EM (T300) 

TASK(TNEM  ,PW  _ _,TRTS)  [user's  name] 

ATTACH , AGO , BINSEMCON , ID=EM71 606 . 

ATTACH ,F1 ,TALWYATT , ID=EM71602 , MR=1 . 

ATTACH , LI BA , AN APAC ,ID=EM71605, MR=1 . 

^ LIBRARY (F1,LIBA) 

MAP (PART) 

LDSET (PRESET=NGINF) 

AGO. 

; TALPLOT  (.TAPED 

i 


i 

I 

{ 

I 

) 

"i 

j 


I' 

I 


t 

I- 
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EXIT. 


7/8/9 

[Data] 

0/6/7/8/9 

If  I PLOT  = 0 and  ITYPE  = 1 on  card  set  1 (plots  come  out  on  the  MDS 
system  using  the  ANAPAC  plotting  software) , then  the  control  cards  for 
running  SEMCON  are  as  follows: 

EM (T300) 

TASK(TNEM ,PW ,TRTS)  [user's  name] 

ATTACH , AGO , BI NSEMCON , ID=EM71606 . 

ATTACH , LIBA , ANAPAC, ID=EM71 605, MR=1 . 

LIBRARY (LIBA) 

MAP (PART) 

LDSET (PRESET=NGINF) 

AGO. 

EXIT.  i 

7/8/9 

[Data]  I 

0/6/7/8/9  1 


10.  CONCLUSIONS 

The  SEMCON  code  has  provided  the  means  to  einalyze  various 
semiconductor  devices  and,  hence,  assist  in  making  EMP  vulnerability 
assessments.  The  code  has  accomplished  this  objective  by  utilizing  very 
straightforward  techniques  to  calculate  the  least-squares  fits  for 
damage  curves  of  semiconductor  devices.  For  EMP  efforts,  the  fits 
generated  by  SEMCON  are  considered  to  be  reliable  when  good,  accurate 
data  are  used.  Also,  the  code  has  proven  to  be  quite  effective  in 
tabulating  and  reducing  numerous  amounts  of  data.  Thus,  SEMCON  is 
considered  to  be  a useful  code  in  semiconductor  damage  analysis  and  is 
written  in  such  a way  as  to  make  it  readily  adaptctble  to  numerous 
computer  systems. 
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ABBREVIATIONS 

Energy 

Failure 

Current  gain  after  device  was  pulsed 
Current  gain  before  device  was  pulsed 
Peak  current 
Device  number 
Power 

Polarity  of  applied  pulse  (i.e.,  reverse  or  forward 
direction) 

Pulse  number 

Failure  time 

Pulse  duration 
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Peak  voltage  | 
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Voltage  after  forward  breakdown  voltage 
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APPENDIX  A. “SAMPLE  PRODUCTION  RUN  OF  SEMCON 


This  appendix  contains  examples  of  plots  availedjle  through  SEMCON 
and  the  computer  center  at  the  Mobility  Equipment  Research  and 
Development  Center,  Fort  Belvoir,  VA. 
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APPENDIX  B. — SAMPLE  PRODUCTION  RUN  OF  SEMCON 

This  appendix  contains  examples  of  plots  available  through  SEMCON 
and  the  Mohawk  Data  Sciences  computing  system. 


pmcsDim  PAOS  blank.not  riiKSD 
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